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Abstract. We show theoretically that nonrelativistic nearly- free electrons in solids 
should experience a trembling motion (Zitterbewegung, ZB) in absence of external 
fields, similarly to relativistic electrons in vacuum. The Zitterbewegung is directly 
related to the influence of periodic potential on the free electron motion. The frequency 
of ZB is uj sa Eg/%, where E g is the energy gap. The amplitude of ZB is determined by 
the strength of periodic potential and the lattice period and it can be of the order of 
nanometers. We show that the amplitude of ZB does not depend much on the width 
of the wave packet representing an electron in real space. An analogue of the Foldy- 
Wouthuysen transformation, known from relativistic quantum mechanics, is introduced 
in order to decouple electron states in various bands. We demonstrate that, after the 
bands are decoupled, electrons should be treated as particles of a finite size. In contrast 
to nearly-free electrons we consider a two-band model of tightly bound electrons. We 
show that also in this case the electrons should experience the trembling motion. It is 
concluded that the phenomenon of Zitterbewegung of electrons in crystalline solids is 
a rule rather than an exception. 
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1. Introduction 

Zitterbewegung (the trembling motion) was theoretically devised by Schroedinger [1] 
after Dirac had proposed his equation describing free relativistic electrons in vacuum. 
Schroedinger showed that, due to a non-commutativity of the quantum velocity 
v = dHo/dp with the Dirac Hamiltonian Hp, relativistic electrons experience the 
Zitterbewegung (ZB) even in absence of external fields. The frequency of ZB is about 
uj = 2itlqc 2 /h and its amplitude is about the Compton wavelength A c = h/moc ~ 
3.86 x 10~ 3 A. It was later understood that the phenomenon of ZB is due to an interference 
of electron states with positive electron energies (E > m$c 2 ) and those with negative 
energies (E < m c 2 ), see [21 EJ H|- In other words, the ZB results from the structure 
of the Dirac Hamiltonian, which contains both positive and negative electron energies, 
and it is a purely quantum effect as it goes beyond Newton's first law. 

An important step in the understanding of ZB was made by Foldy and Wouthuysen 
[5], (see also 0,[7]), who showed that in absence of external fields there exists a unitary 
transformation that transforms the Dirac Hamiltonian into a Hamiltonian in which 
positive and negative electron energies are decoupled. While solutions of the Dirac 
equation are four-component functions, the transformed states for the positive energies 
have only two upper non-vanishing components and those for the negative energies 
have only two lower non-vanishing components. Now the above mentioned interference 
between the positive and negative energy states can not occur and there is no ZB. 
Instead, in the new representation the electron is not a point-like particle, but it acquires 
a 'quantum radius' of the size A c . The interpretation of the two pictures is until present 
not quite clear, see [HJ [91 [101 [TTJ, [131 121 EH]- To our knowledge, the ZB for free electrons 
has never been directly observed. However, in the presence of the Coulomb potential 
the ZB is manifested in appearance of the so called Darwin term [21 El E] • 

It was pointed out some time ago that the Zitterbewegung also may occur in non- 
relativistic two-band systems in solids [15]. It was shown that, similarly to the relativistic 
case in vacuum discussed above, the consequence of the ZB is that it is impossible to 
localize the electron better than to a certain finite volume. Recently, an analogy between 
the Dirac description of electrons in vacuum and the coupled-band k ■ p formalism for 
electrons in narrow-gap semiconductors (NGS) and carbon nanotubes (CNT) was used 
to demonstrate that the ZB should occur in these systems [T6l ITT) . It was shown that, 
in agreement with the 'semi- relativistic' analogy jTS], the ZB frequency is always 
uj ~ E g /h, where E g is the energy gap between the conduction and valence bands. 
The amplitude of Zitterbewegung in NGS and CNT was estimated to be Xz = h/m^u, 
where is the effective electron mass and u ~ 10 8 cm/s is the maximum electron 
velocity in the system. The ZB length in NGS and CNT turns out be 10 — 100A, i.e. 
10 4 — 10 5 times larger than in vacuum. A much lower ZB frequency and its much higher 
amplitude, as compared to vacuum, should make the ZB much more readily observable 
in semiconductors. The Zitterbewegung was also recently proposed in two-dimensional 
systems exhibiting spin splitting due to structure and bulk inversion asymmetry [20J, 
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and in 2D graphite [21J. A phenomenon similar to the ZB was proposed for electrons in 
degenerate valence bands in the presence of an external electric field [22J . Very recently, 
a unified description of the Zitterbewegung of electrons in different solid state systems 
was attempted [23J. 

In view of this recently published work we want to investigate the question of 
whether the phenomenon of Zitterbewegung in solids is a rule rather than an exception 
or vice versa. To this end we consider two limiting models for electrons in solids: 
nearly-free electrons, for which the periodic potential of the lattice may be treated as 
a week perturbation, and tightly-bound electrons, for which the periodic potential may 
not be treated as a perturbation. Since we are interested in general properties of the 
Zitterbewegung, we do not insist on details of the band models in question but rather 
concentrate on essential features that result in this phenomenon. Although we deal with 
non-relativistic electrons in solids, we use methods of relativistic quantum mechanics to 
investigate an alternative picture in which the trembling motion is replaced by a kind of 
electron 'smearing' in real space. The reason, that a somewhat mysterious phenomenon 
of Zitterbewegung of electrons in vacuum has never been observed, seems to be related 
to its very high frequency and very small amplitude. The corresponding phenomenon 
in solids would have much lower frequency and much larger amplitude. The underlying 
hope motivating our work is, that a more thorough theoretical understanding of the 
trembling motion will lead to an experimental detection of the phenomenon. This 
would not only deepen our knowledge of electrons in solids but also represent a great 
success of the relativistic quantum theory. 

Our paper is organized in the following way. In Section II we give the basis of 
nearly-free electron formalism, Section III treats the resulting Zitterbewegung using 
Schroedinger's method of the equation of motion. In Section IV a more realistic 
description of the ZB is presented in which electrons are treated as wave packets. In 
Section V we use the Foldy-Wouthuysen transformation known from the relativistic 
quantum mechanics to obtain an alternative electron picture. Section VI treats the 
Zitterbewegung in case of tightly bound electrons. In Section VII we discuss the obtained 
results and confront them with the previous work. The paper is concluded by a summary. 

2. Nearly- free electrons 

The beginning of this Section is standard, but it is needed for further developments. We 
consider an electron in the presence of an external periodic potential V(r) = V(r + r a ), 
where r a is a translation vector of the lattice. The periodic potential V(r) may be 
expressed by the Fourier series V{r) = J2i Vi exp(iZr), where I are reciprocal lattice 
vectors and Vi are Fourier components of the periodic potential. For a real potential 
there is V,* = V-i- The wave function of an electron has the Bloch form 




1 



(1) 
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where V is the crystal volume and k is the wave vector. Inserting the wave function 
^ k (r) into the Schroedinger equation one obtains the well-known equation for the energy 
E and the coefficients ai 

^-^(fe + 2 )^ = E^- 3) (2) 

where m is the free electron mass and g are reciprocal lattice vectors. 

In absence of the periodic potential there is ai = 0, a = 1, and E = h 2 k 2 /2m = e k 
is the free electron energy. For a weak periodic potential we may treat V(r) as a 
perturbation and approximate ai for I ^ by retaining only linear terms in Vj. We 
obtain then 

ai = (3) 

and a = 1. The perturbed energy is 

E = e k + Y,^^- (4) 

For weak potentials the correction to the free electron energy is small. This, however, 
is true only if e k ^ e fc+i . For k and I = q such that e k = e k+q we expect a q to be 
comparable to a and the potential may not be treated as a weak perturbation. The 
well known way to treat this problem is to use the approximation for nearly degenerate 
levels in which we neglect in fl2]) all ai except ao and a q . We then find 

| (E - e k )a = V- q a q _ 
\ (E - e k+q )a q = V q a . 

Equations (j^D are equivalent to dealing with the Hamiltonian 

which is valid for k such that e k+q ~ e k . Hamiltonian (jH]) has two eigen-energies 

E m = f ^±^ ± yW + ^p} 2 , (7) 
and two eigen-states 
II >= 



where 





E A + A 


i( 


V* 


H 


-v q 




E A + A 


( e fc+<J 


- Cfe) = 


Ml 


2 + A 2 



|2>=-( „ \ ), (86) 

A = \^ +q ~ e k ) = £- ( fcq + ^ ) , (9) 

£a = V / I^l 2 + A2 (10) 
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Figure 1. Energy versus k z for free and nearly- free electrons (schematically). The 
dashed line shows E — h 2 kl/2mo dispersion for free electrons plotted in the first 
Brillouin zone. The solid line indicates E(k z ) dispersion for the free motion weakly 
perturbed by a periodic potential. 



and N = ^2E A (E A + A). 

Taking only one vector q in the reciprocal lattice we are in reality considering a 
one-dimensional problem. To fix our attention we consider the symmetry of a simple 
cubic lattice. As to the two points in the reciprocal lattice, see flHJ), it is convenient 
to take 1 = and I = q = [0,0, — 27r/a], where a is the lattice period. Then the two 
parabolas e k and e k +q in (pE|) cross at k = [0, 0, ir/a], which determines the Brillouin zone 
boundary on the positive k z axis, see figure [D The energy gap at the zone boundary 
is Eg = 2\V q \. The energies A and E A depend in reality only on k z . When using the 
nearly-degenerate perturbation theory based on ([6]) we should keep in mind that this 
procedure is only valid near the degeneracy point k z = n/a, but it progressively ceases 
to work lowered toward zero. 

We note that E A of (TTDT) is analogous to the relativistic dispersion relation E(p) = 
[(m c 2 ) 2 + c 2 p 2 ] 1 ^ 2 . Thus \V q \ 2 = (E g /2) 2 corresponds to (m c 2 ) 2 , while A 2 , which is 
quadratic in momentum, corresponds to c 2 p 2 . 

It is convenient to split the Hamiltonian ([6]) into two parts 

H = H A + H k , (11) 

where 

M£3) 
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and 

H k = \(e k+q + e k ) ^ J J. (13) 

In terms of the Pauli matrices the Hamiltonian ffl2|) reads 

H A = Aa z + M{V q )a x - %{V q )a y . (14) 

Hamiltonian Ha has the form reminiscent of the Dirac Hamiltonian for relativistic 
electrons in vacuum, while the part H k is proportional to the unity matrix and it can 
be treated as a c-number. The decomposition ( Till is directly related to the two terms 
in the energy ([7j). 

The quantum velocity is v = OH / d{hk). We calculate 

v = Va + Vk = u A cr z + u k l, (15) 

where ua = dA/d(hk) = (h/m )q/2 and u k = (h/m )(k + q/2). It follows from (TT5T) 
that the quantum velocity v is an operator, not a number. Since A depends only on k z , 
the only non- vanishing component of Va is v Az - In the following we drop the index z. 

Eigen- values of the quantum velocity va are =F(h/mo)(ir/a). This seems 
paradoxical, as it means that the quantum velocity takes only two constant (and 
extreme) values. A similar result is obtained for the Dirac equation describing relativistic 
electrons in vacuum, for which the eigenvalues of the quantum velocity are ±c. It is 
known that this feature is related to the phenomenon of Zitterbewegung. 



3. Zitterbewegung 

It can be easily verified that the quantum velocity (fl~5|) does not commute with the 
Hamiltonian ([6]). This means that dv/dt does not vanish. The original derivation of 
Schroedinger's is based on the quantum equation of motion (see also [12] ). Let us 
calculate the time dependence of v. We have 
dii 

th— = [v A ,H] + [v k ,H]. (16) 

Since H k and v k are unity matrices, they commute with any number matrices. Therefore 
dv k /dt = l/(ih)[v k ,H] = 0, so that v k (t) = v k0 . Thus 

where the anti-commutator of {va,Ha} = 2uaA. Hence 

. h dvA = 2 y A H A _ 2ma A. (18) 

(JjL 

Let us calculate the second time derivative of va 

ih 3 -^ = i%H A ] = ^[2vaHa ~ 2u A A,H A ] = 
dt l dt in 

= ^-[va,Ha]Ha = 2^Ha- (19) 
in dt 
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This represents a differential equation for dv A /dt. Its solution is 

^ = A eM-2iH A t/h), (20) 

where A is a constant operator. Inserting (1201) to ffT8]) 

iftA exp(-2iH A t/h) = 2v A H A - 2m a A. (21) 
Solving (EH) for -£>a we obtain 

v A (t) = i^i exp(-2zF A t/^)i7X 1 + u A AH A \ (22) 

Integrating ( 1221) with respect to time and adding the term due to z component of Vk we 
finally find 

z(t) = z + v kz0 t + u A AH A ~ 1 t + 

- ^h 2 A (exp(-2iH A t/h) - l) H A 2 . (23) 

In order to find A we use (1211) for t = 

Ao4K"A)^U, o)' (24) 

At t = there is z(0) = z . Similarly, it follows from (1221) and (1241) that 'Da(O) equals to 
the z component of the initial velocity v A from (fTBT) . 

In order to interpret the result (1231) we observe that the eigen-energy of _£/a is 
±£ A , where E A is given by (TTD]). There is if^ 1 = H A /E A 2 and F^ 2 = 1/^a 2 - The 
exponential term in (123]) is (see Appendix) 

f-2iH A t\ (2E A t\ .Ha . (2E A t\ 

exp —fr~ = cos {— ) ~ 1 e: sm InrJ • (25) 



The first three terms in ([231) describe the classical motion. The last term, 
according to (1251) . describes oscillations with the frequency uj = 2E A /h. This frequency 
corresponds directly to the interband energy 2E A , as seen in figured! 

Since the Hamiltonian H A is a matrix, the position z(t) is also a matrix. We have 
explicitly 

z(i) = Uw *.(«)J' ( 1 

where 

. r .s ^a|^| 2 . f2E A t\^ u A A 2 t 

zn(t) = 3 sm 1^^— J + ^ 2 + + zq. (27) 

The component z 2 2 it) has negative signs of the first two terms. Further 



U ; 2E A 



(2E A t\ 



A /2E A t \ 1 u A AO , , 
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where z = z(0) and Suit) = z 2 i(t)*. 

The amplitude of the oscillating term in (T27D is hu^\V q \ 2 /2E A z « nh 2 / (2m a\V q \) = 
Xz/2, where the Zitterbewegung length is defined as 

A Z = — ~~7Trp (29) 
m a\ V q \ 

It corresponds to the Compton wavelength in relativistic quantum mechanics. Its 
numerical estimation is given below. 

In agreement with the history of the subject, as described in the Introduction, we 
can legitimately call the above oscillations the Zitterbewegung. We shall discuss the 
subject of ZB more thoroughly below. Here we emphasize how little we have assumed 
to obtain the trembling motion - we have only perturbed the free electron motion by a 
periodic potential. 

4. Wave packet 

Now we consider a more realistic picture describing an electron in terms of a wave 
packet. Taking, as before, q = [0,0, — 2n/a] and using the fact that A and depend 
only on k z , we take a gaussian packet in the form 

^(*) = ^r"i;a / ex P ( -\ d2 { k z - k z0 ) 



2 

exp(ik z z)dk z ^ q J > ( 30 ) 

where d determines the packet's width and k Z Q fixes its center in the k z space. Averaging 
the oscillating part of the motion over the wave packet we obtain from (1261) 

<^(z)\z° sc (t)mz) >= 

: J zl?{t) exp (-d 2 (k z - k z0 ) 2 ) dk z , (31) 




where z°f c (t) is given in 

To begin, let us take the packet to be a delta function in the k z space centered 
at k Z Q. This corresponds to a completely non-localized packet in the real space. We 
can then take various k Z Q values beginning with k z o = i\ja at the zone boundary. In 
figure [2] we show the calculated ZB oscillations of zn(t) for different values of k z o- It 
can be seen that, as k z0 diminishes from n/a toward the zone center, the amplitude of 
ZB quickly drops. This should not be surprising since, as is well known (see figured]), 
the effect of the periodic potential on the free electron motion is the strongest at the 
zone boundary k z = ir/a, where the minimum gap occurs. (We do not consider here the 
gap at k z = for the upper branch). Figure [3] shows the amplitude of ZB for d — oo 
and d = SOr^ 1 , as calculated from fl3T|) . The quantity r B = 0.53A is the Bohr radius. 
For d = oo we deal in ( 1301 and ( 13TT) with a delta-function and the solid line in figure [3] 
follows the dependence Uu&\V q \ 2 / (2E& 3 ) of (|27|) . When the width of the wave packet 
increases (d decreases) the amplitude of ZB for k Z Q w rr/a diminishes and, as k z o is 




Figure 2. Zitterbewegung oscillations of nearly-free electrons versus time, calculated 
for a very narrow wave packet centered at various k z o values. The band parameters 
correspond to GaAs, see text. 
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Figure 3. Amplitude of Zitterbewegung of nearly- free electrons versus packet center 
k z0 , calculated for two widths of the wave packet. Symbol denotes the Bohr radius. 
Material parameters correspond to GaAs. 
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Figure 4. Zitterbewegung contribution to the velocity for nearly-free electrons versus 
time, calculated for a very narrow wave packet centered at various k Z Q values. The 
band parameters correspond to GaAs, see text. 



lowered, it becomes independent of the width. Since our model is not valid for k Z Q near 
zero, we are limited in our considerations to not too small d and h z o values. 

However, since the ZB amplitude diminishes so quickly with diminishing k Z Q, it is 
justified to limit the considerations of ZB to the vicinity of the band extremes. It can 
be seen from figure [2] that, with decreasing k z o, the frequency of ZB increases. This 
increase follows from the behavior of the gap 2E A , as illustrated in figured! 

In order to calculate numerical values of the ZB we need to specify material 
parameters. As a matter of example we take V q = E g /2 = 0.76eV and a = 5.6A, 
corresponding to GaAs. This gives Xz = 5.6A. This value can be compared with 
Xz = 10 — 13A for GaAs, as obtained with the use of kp theory for the fundamental gap 
in GaAs at k = [16]. The above estimation of Xz based on the simple model is better 
than one could expect. Clearly, if we take the V q value corresponding to E g = 0.23eV 
for InSb, Xz would be seven times larger. 

Next we calculate an observable quantity that is the electric current caused by ZB. 
It is given by the velocity multiplied by the charge. The oscillatory part of velocity is 
given by the first term in (1221) . We average it using the wave function (1301) which selects 
the component Vn(t) 

. M u A \V q \ 2 {2E A t\ u A A 2 

M*) = -^cos(- r j+-g r . (32) 

The results for the velocity, computed with the gaussian wave packet (1301) . are plotted 
in figure HI They are quantitatively similar to those shown in figure [2], the ZB frequency 
is clearly the same, the amplitude decreases with increasing k Z Q. The phase, however, is 
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different and the velocity v\\(t) is not zero at t = 0, but is equal to «a from (IT51) . Also, 
the velocity does not oscillate around zero, it oscillates around the value resulting 
from the second term in ( |32i) . For k Z Q = it /a we can find analytically 

„g) = UA {1 - v^Cexp(C 2 )[l - erf(C)]}, (33) 

where ( = d/Xz and erf(x) is the error function. For narrow packets (large d) the shift 
tends to zero. 



5. Foldy-Wouthuysen transformation 

As mentioned in the Introduction, Foldy and Wouthuysen [5] proposed a transformation 
which, in absence of external fields, transforms the Dirac Hamiltonian for relativistic 
electrons in vacuum into a form in which positive and negative electron energies are 
separated. It was recently shown by Zawadzki [161 02] that similar transformations exist 
for the A; • p Hamiltonians describing band structures in narrow gap semiconductors and 
carbon nanotubes. Since also the Hamiltonian H A of (jT2l) bears similarity to the Dirac 
Hamiltonian, we can expect that a similar transformation exists for nearly-free electrons 
as well. This is indeed the case. 

We define a unitary transformation 

E A + (3H A 




/2E A (E A + A) 

where $ = ( " I . It is easy to verify that UW = 1. Further 



(34) 



UH A U ] = E A (3 (35) 

and obviously UHfJj^ = Hf~ since Hk is proportional to the unity matrix (see ffTB"]) ). 
Thus, for the transformed Hamiltonian the eigen-energy problem factorizes into two 
independent problems for positive and negative E A energies. This means that the wave 
function corresponding to the positive energy has the lower component equal to zero 
while the wave function for the negative energy has the upper component equal to zero. 
However, this is true also for other wave functions in the transformed representation, as 
we show below. 

We consider an arbitrary wave function *&(z) in the two-component representation. 
It can be expressed in general in the form 

V(z) = / u(k z ) exp(ik z z)dk z = + (36) 

where 

*±W = \ r f 1 ± «(*!) exv(ik' z z)dk' z . (37) 



2 J— oo V E A 
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Let us now transform the above functions using the U operator: *ff'±(z) = U^±(z). 
After some manipulations we obtain 

1±$ f 00 ( E A xV2 



x (l±^\u(k' z )exp(ik' z z')dk' z . (38) 

The /^-dependent function under the square root may be put under the integral sign, 
as explained in [5]. 

The functions ^f' ± (z) have the above mentioned property of having only the upper 
or the lower non- vanishing components, which is guaranteed by the pre-factors (1 ± p). 
Using the inverse Fourier transform 



u ( k ' z ) = r ^{z')exp{-ik' g z')dz, (39) 

/oo 
K^{z,z)^{z)dz\ (40) 
-oo 



\ 1/2 



~±, >, l±P 1 f°° ( E A 



we have 



where 



2 2n J-oo \2(E a + A) 
(l ± ?±\ ex P «(z - (41) 

The kernels /C ± (z, z ) are not point transformations. To illustrate this we will transform 
the eigen-function of the position operator z in the old representation, i.e. the Dirac 
delta function multiplied by a unit vector 

*(z) = 5(z - z ) ( M. (42) 







The transformed functions are 



and 



where 



and 



<(z) = /C+(z,z )( I J (43) 



9L(z)=K^(z,zo)(° 1 ), (44) 



^(^o)^£/l + A e <(-)^ (45) 



z ) = —% \ = dk z . (46) 

27TV2 J -co JE A (E A + A) 
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It is seen explicitly that the transformed functions *ff + (z) and ^'_(z) have vanishing 
lower or upper components, respectively. Both K.xi(z,Zq) and K.2i(z,z ) are normalized 
to delta functions. To demonstrate that the transformed functions are characterized by 
a certain width, we calculate their second moments A4f. Since fCf^z, zq) = £u(C); 
where ( = z — zq, we have 



(££(C))Wi(C)dC 



AC 



\ 



1 + 



A 



k'. 



Ea , 



x 



\ 



A,. 



Ea,, 



C ! 2 dC ! dk z dk z 



(47) 



In (I47p we used subscripts k' z and k" z to indicate the variables of the integration. Using 
relations Q = i(d/dk' z ) exp(—ik' z () and similarly for k" z the triple integral is reduced to 

/ 1 



1 

2^ J- 



d 
dk' 



\ 



\ 



1 + 



A 



dk. 



32' 



(48) 



where Xz is defined in (1291) . The above result has been calculated directly using A and 
from (Q and ffTDl . Note that since ^^(C) has the dimension [m _1 ], the dimension 
of M2 is [m]. 

Similar calculations for the second moment of K,2i{0 gi ye -^2 = A^/32. Thus 
the transformed functions for the upper and lower energies are characterized by the 
same widths, as should be expected. We will discuss physical implications of the above 
calculations in Section VII. 

When transforming various wave functions from the two-component representation 
to the one-component representation with the use of (|40l) . it is important to know more 
about the kernels lC ± {z — z), as given by fT4Tl) . As an example we will calculate and 
plot /C^(C) given by (gg) 

Because both A and Ea are centered at k' z = —q z /2 = (tt/a), it is convenient to 
change the variables k' z — » k z — q z /2. Then we obtain 



/C+ (C) = expHg,C/2) x K+((), 



where 



^i + i(0 




exp(ik z Qdk z 



(49) 



(50) 



and _Ea C 



A c 



2m 



k z q z 



(51) 



V q \ 2 + (A ) 2 . The quantities A and Ea are centered at k z = 0. After 



the change of variables, we singled out the rapidly oscillating part of /Ci~i(C)> which is 
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related to the position of bands extremes in the k space. The remaining part i^n(C) is 
a smoothly varying function of ( with a singularity at ( = 0. 
We consider the integrand in ( 1501) 




(52) 



For k z — > oo the function B(k z ) tends to unity as 1 — 0(k^), while for k z — > — oo it tends 
to zero as 0(k z ). Therefore the integral ( 1521) is poorly convergent. Nevertheless, we can 
calculate it with the help of the Heaviside function Q(k z ), which has a similar behavior 
to the integrand ( 1521) for k z — > ±oo. We have 

1 / f 00 

Z7T \J-oo 

+ Q(k g )exp(ik g Qdk e ) . (53) 

J — oo / 

The second integral is /f^ 0(fc z ) exp(ik z ()dk z = i/( + 7r8((). Thus the real and 
imaginary parts of ifn(C) are 

^i + i(C)] = r- / - (^)] 008(^0*- + o 5 (0 (54) 

Z7T J-oo Z 

and 

1 r°° 1 

S[^i + i(C)] = 5- / " (^)] sin(feC)dfc- + Try (55) 

Z7T J-oo Z7TC, 

The above integrals are carried out numerically. The results are plotted in figure [51 Both 
5i[-Kii(C)] an d 9f[ifji(C)] are singular at ( = and for large £ they decay exponentially. 

Since we are dealing with weakly convergent integrals it is of interest to verify the 
accuracy of the above numerical calculations. To this end, we check two sum rules 
holding for defined in ( 1501) (see [2]). The first rule is 

1 / A°~~ 1 







v/2' (56) 

Now we verify that the same result is obtained from the numerical calculations of 
integrals in (154j) and (l5"5l) . Since the imaginary part of i^n(C) is an odd function of 
C, it gives no contribution to Sq. The numerical calculation of the real part of 
gives l/v2 with the accuracy of 1CT 3 . 

The second sum rule is related to the second moment of Calculations 
similar to the presented above give 




(57) 



J — oo 

_L ^ 

Taking into account the normalization (|56|) we obtain the extension of function to 
be A^/2, which is in exact analogy to the relativistic Dirac electrons [2], see also [16J. 
Similar results are obtained for the Sq and S% integrals defined using K%\ function. 



Zitterbewegung of electrons in semiconductors 



15 




Figure 5. Real and imaginary parts of the kernel of Foldy-Wouthuysen 
transformation for nearly- free electrons (the smooth component) versus £ = z — zq. 
Material parameters correspond to GaAs. 



The FW transformation is not the only transformation which can decouple positive 
and negative energies in the field-free case. In the relativistic quantum mechanics 
described by the Dirac equation other transformations were devised, see for example 
Cini-Touschek [24]. In a recent paper Mulligan [25] introduced a still different 
transformation separating the 4x4 Dirac equation into two 2x2 equations for the 
electron and antielectron, respectively. It is possible that an analogous transformation 
would be possible for the nearly-free non-relativistic electrons considered above. 

6. Tightly bound electrons 

In the preceding sections we considered the case of a weak periodic potential acting on 
free electrons and we showed that this potential leads to the ZB. Now we consider an 
opposite limit of strong periodic potential. 

An effective treatment of a strong periodic potential is the tight-binding method. 
We use here as an example the so called Empirical Tight Binding Method. In this 
model one takes one s orbital per cation and three p orbitals per anion including nearest- 
neighbor and second nearest-neighbor interactions. Spherical approximation is assumed, 
so all sp and p bands are isotropic and their fc-dependence is given by the F" — X 
dispersion. The model was used for calculating magnetic interactions in dilute magnetic 
semiconductors, approximating the band structure of Cdi-^Mn^Te within the whole 
Brillouin zone [26]. This scheme provides a good semi-quantitative description of both 
upper valence bands as well as the lowest conduction band. In the basis x a ,y a , z a , s c , 
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(c-cation, a-anion) the Hamiltonian is 



H 



k 




iV k 





*(3) 
l k 









,(3) 



iV k \ 





t 



(1) 



(5f 



where 



a(i) 



,(2) 



»(3) 



V, 



e c + 4C[l + 2cos(-aA;)], 



= e a + 4A 2 + 8A cos(-a/c)], 

1 1 

= e a + 4Ai[l + cos(-afc)] + 4A 2 cos(-aA;), 

AVcaSmi-ak). 



(59) 
(60) 
(61) 
(62) 

Six parameters in (J59>(|g2]): e c = 3.16eV, e a = O.leV, V ca = 1.103eV, C = 0.015eV, 
A\ = 0.13eV, A 2 = 0.15eV are the Slater-Koster parameters in notation used in |26j, 
and a = 6.482A is CdTe lattice constant. The parameters e c (e ) are cation (anion) 
on-site energies, V ca is a single nearest- neighbor hopping parameter, while C, A%, A 2 are 
the second-neighbor parameters. 

Hamiltonian (1581) can be factorized giving one doubly degenerate energy band 
E 2 ,3 = ts(k) coming from p y and p z orbitals, and two energy bands coming from the 
s — p x interaction. The sp x Hamiltonian can be written in the form 



H. 



sp 






(63) 



>\"k +^r ; )- The Hamiltonian (1631) is very similar to 



where A k = |(4 } — 4 ) and r fc = §(*. 
that for nearly- free electrons of (jTTj) - (jl~3"]) . The difference is, that A k and V k depend on 
the absolute value of |fc|. Also, now the minimum band gap occurs at k = 0. All eigen- 
energies of the Hamiltonian H sp have periodicity of the lattice constant. The quantum 
velocity v = dH sp /(dhk) does not commute with H sp . To calculate f(t) we use the 
Heisenberg picture 

f{t) = ex-p(iH sp t/h)r(0) exp(-iH sp t/h), (64) 

which gives for the f(t) matrix 



fn(t) 



+ 



'V k A' k V k -V k A k V k ' 

2E A 3 
' A k V k Vl + AjA'A 



sm 



{2E A t\ 



h 

v v t + r . 



) + 



(65) 
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i 




k Q [n/a] 

Figure 6. Amplitude of Zitterbewegung of tightly bound electrons versus packet 
center fco, calculated for three widths of the wave packet. Symbol tb denotes the Bohr 
radius. Material parameters correspond to CdTe, see text. 



The component f 2 2{t) is given by (165]) with the changed signs of the first two terms. 
Further 



A , , A k V k - A k V k 

hi(t) = 



_k 1 

2E, 2 



{2E A t\ 



A 



and r i2 = r 21 . In ( |65|) - (!66|) the prime denotes a differentiation with respect to k, 
and t>r — dT k /(dhk). Equations fl65l) - (|66l) are formally similar to (!27l) - (l28l) . the main 
difference is that now there appear terms related to the dependence of V k on k. 

We calculated the ZB oscillations of fn(i), as given by fl65l) . using the gaussian 
packet of the form f[5U|) . The results are similar to those illustrated in figure [2J In figure 
[6] we show the calculated amplitudes of ZB as functions of the packet center k for three 
widths of the packet. Here the fc -dependence does not have a maximum at k = 0, 
because at this point the interaction V k between the bands vanishes, see fl62l) . 

The presence of a strong periodic potential leads to two effects. First, the quadratic 
dispersion relation oc k 2 for free electrons is replaced by a periodic one. Second, the 
potential mixes s and p x states to form two sp x energy bands. This mixing leads to the 
Zitterbewegung. 
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7. Discussion 

The main result of our work is that both in the case of nearly-free electrons and in the 
opposite case of tightly bound electrons we predict the Zitterbewegung phenomenon. 
Comparing this result with the previous work, in which the ZB was predicted with the 
use of LCAO [15] and k p theory [161 H7j, we conclude that ZB is not due to a particular 
approach. In fact, the mathematics is quite similar in all above theories and, although 
we deal with non- relativistic electrons, it resembles the formulation of relativistic 
quantum mechanics for free electrons in vacuum. It is clear that the fundamental 
underlying reason for the appearance of ZB in solids is the periodic potential of the 
lattice. Particularly instructive in this respect are figures |2] and [3] of the present paper 
which show that the amplitude of ZB is directly related to the effect of periodic potential 
on the free electron motion. 

This result is not surprising. Without specifying any particular band model we deal 
in solids with the periodic Hamiltonian H = p 2 /2niQ + V(r) and the electron velocity 
is v = dH /dp = p/m . It follows that dv/dt = l/(ih)[v, H] ^ 0, i.e. the velocity is not 
constant in time because of the periodic potential V(r). In this perspective the various 
models mentioned above simply illustrate how this result comes about. 

Clearly, the velocity does not commute with the Hamiltonian in the presence of 
other potentials as well. However the periodic potential is special because, due to the 
Bloch theorem, the electrons can propagate in the perfect crystal without scattering 
and the quasi-momentum hk is a good quantum number. For this reason it is possible 
to treat the electrons as almost free particles and replace the influence of the periodic 
potential by an effective electron mass. Still, as demonstrated in Refs. [161 HZ] an d 
the present paper, within two-band models the basic non-commutativity of v and H 
mentioned above remains in the form of non-commuting 2x2 matrices with the resulting 
Zitterbewegung. If the bands are completely separated, we have H e ff = p 2 /2m* and 
v = p/m*, so that H e ff and v commute and the ZB disappears. However, there is a 
price to pay for this separation. It is shown in Refs. [HH [TTJ and in the present paper 
that, once the electrons are described by a one-band equation (so that their energy is 
completely specified), they should be treated as objects of a finite size. The last effect 
is observable in the presence of an external potential due to appearance of the so called 
Darwin term for free relativistic electrons [2j [3j HJ CEO] as well as semiconductor electrons 

As emphasized throughout our paper (see also Refs. [HUE]), in the two-band model 
the ZB of non-relativistic electrons in solids is in close analogy to the ZB of relativistic 
electrons in vacuum, as first proposed by Schroedinger. The Hamiltonians for the two 
cases are very similar and in both systems the ZB results from an interference of electron 
states corresponding to positive and negative electron energies (2j[3jll]. If, with the use 
of Foldy-Wouthuysen transformation, the states of positive and negative energies are 
separated in the Hamiltonian and in the wave functions, the ZB does not occur because 
the positive energy state (say) has nothing to interfere with. This corresponds to the 



Zitterbewegung of electrons in semiconductors 



19 



separation of bands mentioned above and the conclusions of the two reasonings agree. 

Thus we are confronted with the following choice: 1) We use a two- band description, 
the electrons are point-like particles and they experience the Zitterbewegung. 2) We use 
a one-band description, the electrons do not experience the Zitterbewegung but they 
are characterized by a quantum radius of the size equal to the ZB amplitude. 

The last point is illustrated by fl57|) describing the average 'smearing' of the 
transformed delta function. It is equal to the amplitude of ZB given by ff271) . One 
can say that the separation of energy bands by the Foldy-Wouthuysen transformation 
is equivalent to a certain averaging of the ZB motion. 

It was observed [IT] that, since the Zitterbewegung had been predicted for plane 
Dirac waves, it is not quite clear what the trembling motion means for an electron 
uniformly distributed in space. In this connection it is important that the amplitude 
does not vary much when the electron is represented by a wave packet localized in real 
space, see figure |3j 

Passing to more specific points of our treatment, we emphasize again how little we 
had to assume to derive the ZB in case of nearly-free electrons - it was enough to perturb 
the free electron motion by a periodic potential. This case has certain particularities. 
In the 'typical' two-band situations, both in vacuum J2J HI [10] and in solids [TH ITT] , 
there exists a maximum velocity in the system which plays an important role in the 
theory. In case of nearly-free electrons there is no maximum velocity since the perturbed 
energy branches tend asymptotically to the free electron parabola E = h 2 k 2 /2mo, see 
figure UJ This is reflected in the velocity i)k, while the velocity has the typical 
'two-band' behavior and it is responsible for the ZB. In fact, the velocity has the 
maximum value. It is equal to (h/ mo) (tt /a), see our eigen- value considerations after 
ffl5]) . For a = 5.6A the maximum velocity is 6.4xl0 7 cm/s, which should be compared 
with u = 1.3 x 10 8 cm/s obtained from the k ■ p theory for GaAs and other III-V 
compounds [16j. Again, our simple model gives quite a reasonable estimation. In the 
nearly-free electrons model Xz is proportional to l/V q ~ 1/E g (see (1291). which agrees 
with the k ■ p approach |16j . where Xz ~ l/ m o ~ V-^s- m narrow gap materials Xz can 
be as large as tens of angstroms. 

As far as the phenomenon of ZB is concerned, the behavior of nearly-free and 
tightly bound electrons is quite similar. The main difference comes from the fact that 
in nearly-free case the Fourier coefficients of the periodic potential do not depend on 
the wave vector k, whereas in the tightly bound case they go to zero for vanishing 
k. As a result, in the first case (see figure [3]) the ZB amplitude is highest for the k zQ 
corresponding to the minimum energy gap, while in the second case (see figure El) the 
maximum amplitude is shifted with respect to the minimum gap. For the tightly bound 
case the electrons in t^ bands (see flKTJ ) would not exhibit the ZB, but we do not insist 
on this point since these two bands are not realistically described by the model. 

In order to illustrate that, after the Foldy-Wouthuysen transformation is carried 
out, the former Dirac delta function is 'smeared' into the kernel /C^ of ( 1451) . we not 
only calculate it numerically (see figure [5]), but also calculate its second moment M.\ 
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(see (jUD) and the sum rule S£ of (1571) . The second moment has the advantage of 
using the standard quantum mechanical probability distribution (K^YK^, see ( HTj) . 
Its disadvantage is that this probability distribution is normalized to the Dirac delta 
function and not to a number. As to the sum rules using the single function, their 
advantage is that both the normalization ( 1561) and the sum rule S£ ( 1571) are numbers, 
and (1571) can be simply interpreted as a square of 'smearing' (This procedure was used 
for free relativistic electrons by Rose [2]). The disadvantage of using the single is 
that we have to separate out its 'smooth' part, see ( )49|) . The oscillatory part exp(— iq z C) 
would not appear if we considered the nearly- free electron gap at k z = 0. 

Methods to observe the Zitterbewegung should clearly be adjusted to the 
investigated materials, but it seems that an appropriate tool would be the scanning 
probe microscopy which can produce images of coherent electron flow [271 EE]- This 
technique uses a sharp mobile tip which can sense the electron charge. If one used 
dilute magnetic semiconductors of CdMnTe type [26] , one could employ magnetic effects 
caused by electron oscillations. 

The second category of possible observable effects is related to the problem of what 
happens when electrons are confined to dimensions smaller than Xz- In relativistic 
quantum mechanics one finds statements that a measurement of the position of a 
particle, if carried out with greater precision than the Compton wavelength, would 
lead to pair production [8J. It is clear, however, that the pairs created this way can 
only be virtual, otherwise their recombination would lead to the production of energy 
out of nothing. The virtual carriers could be observed in screening or in magnetism. 
An alternative point of view states that a stiff confinement is equivalent to an infinite 
potential well having the width Ar < Xz, and the electrons will simply occupy the 
lowest energy level in such a well. It is, however, certain that if an electron is confined 
to dimensions Ar < Xz, its energy (or its uncertainty) is of the order of the gap E g 
between the positive and negative electron energies, which means that the one-band 
description is not adequate. We are then back to the two-band model, which was our 
starting point. 

There remain many unanswered questions concerning the trembling motion but it 
appears that in crystalline solids it represents a rule rather than an exception. According 
to the theory, the ZB in semiconductors has decisive advantages over the corresponding 
effect in vacuum. Thus an experimental detection of the trembling motion in solids 
appears to be a matter of near future. 

8. Summary 

We considered theoretically non-relativistic nearly-free electrons in solids for which the 
periodic potential of the lattice may be treated as a weak perturbation on the free 
electron motion. Using the two-band model we showed that electrons experience the 
trembling motion (Zitterbewegung) in absence of external fields, similar to that for free 
relativistic electrons in vacuum. The frequency of ZB and its amplitude were derived. 
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The frequency is u ~ E g /h where E g is the energy gap between the two bands. The 
amplitude Xz depends on the strength of periodic potential and the lattice period. For 
typical parameters Xz can be of the order of 10A to 100A that is 10 4 - 10 5 times 
larger than in vacuum. The trembling motion is also considered for nearly-free electrons 
represented by wave packets, it is shown that the amplitude is not strongly dependent 
on packet's width. The Foldy-Wouthuysen type of unitary transformation, known from 
relativistic quantum mechanics, is used to separate the energy bands. Consequences 
of the FW transformation are investigated. It is demonstrated that, if one uses a one- 
band description, the electrons do not experience the trembling motion but they should 
be treated as particles having size Xz- Tightly bound electrons are considered as well 
to provide an opposite case to nearly-free electrons. Within the two-band model the 
trembling motion is obtained also in this case demonstrating that the ZB phenomenon is 
not related to a specific theoretical approach. It is concluded that the trembling motion 
is directly related to the effect of the periodic potential on the electron and, as such, it 
should occur in many situations in solids. 
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Appendix A. 

First, we prove some properties of H A of ( fl2l) . The eigen- values of H A are ±E A . If 
|1 > and |2 > are eigen-states of H A , then A = \i >< i\, (i = 1,2) are two projection 
operators. It follows that A + A = 1 and H A = E A {P\ — A)- Accordingly 



(A - A) = 




(A.l) 



Since H A l = (1/£ A )(A- A), we have if A 1 = H A /E A 2 . Then H^ 2 = H^H^ = 1/E A 2 
because H\ = E A 2 . For a real b there is 

exp(iH A b) = exp(iE A b)P 1 + exp(-iE A b)P 2 . (A.2) 

Because exp(±iE A b) = cos(E A b) ± ism(E A b) we have 

exp(iH A b) = [cos(E A b) +isin(E A b)}P 1 + 

+ [cos{E A b) - % sin{E A b)]P 2 . (A.3) 

Grouping the terms with cosine and sine functions we get 

exp(iH A b)= cos(E A b)(P 1 + P 2 )+ism(E A b)(P 1 -P 2 ) = 

fi 

= cos{E A b) + i— - sm(E A b). (A.4) 
E A 

This identity is used in ( ]2"5|) . 
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